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O ' Abstract 

o : 

fvq , Wireless network topologies change over time and maintaining routes requires frequent updates. Updates are 

costly in terms of consuming throughput available for data transmission, which is precious in wireless networks. 
\ In this paper, we ask whether there exist low-overhead schemes that produce low-stretch routes. This is studied by 

■ using the underlying geometric properties of the connectivity graph in wireless networks. 
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I. Introduction 

A major challenge in the design of wireless ad hoc networks is the need for distributed routing algorithms that 
consume a minimal amount of network resources. This is particularly important in dynamic networks, where the 
topology can change over time, and therefore routing tables must be updated frequently. Such updates incur control 
traffic, which consumes bandwidth and power. It is natural to ask whether there exist low-overhead schemes for 
dynamic wireless networks that could produce and maintain efficient routes. In this paper we consider dynamically 
changing connectivity graphs that arise in wireless networks. Our performance metric for the algorithms is the 
average signaling overhead incurred over a long time-scale when the topology changes continuously. We design a 
routing algorithm which can cope with such variations in topology. We maintain efficient routes from any source 
to any destination node, for each instantiation.!] of the connectivity graph. By efficient, we mean that we want to 
guarantee that the route is within a (small) constant factor, called stretch of the shortest path length. In order to 
route to a destination, we need only the identity of the destination and not its address i.e., the control traffic to 
maintain the mapping between node identity and address/location is incorporated into the overhead. Therefore, in 
the wireless routing terminology, we have included the "location service" in the control signaling requirement, and 
therefore hope to characterize the complete overhead needed to maintain efficient routes. 

In order to develop and analyze the routing algorithms we utilize the underlying geometric properties of the 
connectivity graphs which arise in wireless networks. This geometric property is captured by the doubling dimension 
of the connectivity graph. A graph induces a metric space by considering the shortest path distance between nodes 
as the metric distance. The doubling dimension of a metric space is the number of balls of radius R needed to cover 
a ball of radius 2R. For example a Euclidean space has a low doubling dimension as will be illustrated in Section 
HH A metric space having a low (constant independent of the cardinality of the metric space) doubling dimension is 
called "doubling". We show that several wireless network graphs (under conditions given in Section [TT]) are doubling 
and therefore enable the design and analysis of hierarchical routing strategies. In particular, it is not necessary to 
have uniformly distributed nodes with geometric connectivity for the doubling property to hold, as illustrated in 
Figure [2] in Section ITT] Therefore, the doubling property has the potential to enable us to design and analyze 
algorithms for a general class of wireless networks. Moreover, for a large class of mobility models, the sequence 
of graphs arising due to topology changes are all doubling (for specific wireless network models). Since there are 
only "local" connectivity changes due to mobility, there is a smooth transition between these doubling graphs. We 
can utilize the locality of topology changes to develop lazy updates methods to reduce signaling overhead. 

We show that several important wireless network models produce connectivity graphs that are doubling. In 
particular, we show that the geometric random graph with connectivity radius growing as ylogn with network size 
n; the fully connected regime of the dense or extended wireless network with signal-to-interference-plus-noise ratio 
(SINR) threshold connectivity; some examples of networks with obstacles and non-homogeneous node distribution. 
We define a sequence of wireless connectivity graphs to be smooth if each of the graphs is doubling and the shortest 
path distance between two nodes in the graph changes smoothly (defined in Section JT]). These for mild regularity 
conditions on the mobility model. 

Our main results in this paper are the following, (i) For smooth geometric sequence of connectivity graphs, we 
develop a routing strategy based on a hierarchical set of beacons with scoped flooding. We also maintain cluster 
membership for these beacons in a lazy manner adapted to the mobility model and doubling dimension, (ii) We 
develop a worst-case analysis of the routing algorithm in terms of total routing overhead and route quality (stretch). 
We show that we can maintain constant stretch routes while having an average network-wide traffic overhead of 
0(n log 2 n) bits per mobility time step. The load-balanced algorithm would require 0(log 3 n) bits per node, per 
mobility time. Through numerics we show that the theoretically obtained worst-case constants are conservative. 

A. Related Work 

Routing in wireless networks has been a rich area of enquiry over the past decade or more. The two main 
paradigms for routing have been geographic routing and topology based routing. Geographic routing (see for 
instance [KKOO] and references therein) exploits the inherent geometry of wireless networks, and bases routing 

'We assume inherently that the round-trip time (RTT) of a packet from source to destination is much smaller than the time-scale of 
topology change. 
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decisions directly on the Euclidean coordinates of nodes. Their performance depends on how well the Euclidean 
coordinate system captures the actual connectivity graph, and these approaches can therefore fail in the presence of 
node or channel inhomogeneity (like in Figure |2] in Section JT]). Another important, but often overlooked, issue with 
geo-routing is that geographical positions of the nodes need to be stored and continuously updated in a distributed 
database in the network, to allow sources of messages to determine the current position of the destination. This 
database is called a location service (see for instance [LJDC + 00]) and must be regularly updated so that source 
nodes can query it. Location services typically rely on some a-priori knowledge of the geographical boundaries of 
the network. This is necessary because these approaches typically establish a correspondence (for example, through 
a hash function) between a node identifier and one or several geographical locations where location information 
about that node is maintained. An important feature of our work is that we consider the total overhead incurred by 
the update and lookup operations of the location service, and the overhead of the routing algorithm itself. 

Topology based routing schemes (see [PR97] and [JMB01]) do not utilize the underlying geometry of wireless 
connectivity graphs, but instead compute routes based directly on that graph. To reduce overhead, most of these 
schemes only establish routes on demand through a route discovery operation, rather than continuously maintaining 
a route between every pair of destinations; in this respect, they differ significantly from their counterparts for the 
wired Internet (such as OSPF, IS-IS, and RIP). Recently established routes are cached in order to allow their reuse 
by future messages. In distance-vector based approaches (e.g., [PR97], this cached state resides in the intermediate 
nodes that are part of a route, whereas in source -routing approaches (e.g., [JMB01]), the cached state resides 
in the source of a route. Despite such optimizations, topology-based approaches suffer from the large overhead 
of frequent route discovery operations in large and dynamic networks. This issue was, in fact, the reason why 
geo-routing approaches have reached prominence. 

Two schemes that utilize the underlying geometry of graphs in static wireless networks algorithms are the works 
presented in [RRP + 03] and the beacon vector routing (BVR) introduced in [FRZ + 05]. Both these schemes are 
heuristics which build a virtual coordinate system over which routing takes place. They were shown to work 
well through numerics. However, they utilize an external addressing scheme to make a correspondence between 
addresses and names. In [TDGW07], routing on dynamic networks using a virtual coordinate system was studied. 
For large scale dynamic wireless networks, these heuristics pointed to significant advantages to using some geometric 
properties for routing and addressing. These results motivated the questions studied in this paper. 

There has been a vast amount of theoretical research on efficient routing schemes in wired (i.e., static) networks 
(see for example [GavOl]). Most of this work has been focused on the trade-off of memory (routing table size) and 
routing stretch. There are two main variants of such routing schemes (i) labeled (or addressed) routing schemes, 
where the nodes can be assigned addresses so as to reflect topological information; (ii) named routing, where nodes 
have arbitrary names, and as part of the routing, the location (or address) of the destination needs to be obtained 
(similar to a location service). This examines the important question of how the node addresses need to be published 
in the network. Routing in graphs with finite doubling dimension has been of recent interest (see [KRX06], and 
references therein). In particular [Tal04] showed that one could get constant stretch routing with small routing table 
sizes for doubling metric spaces, when we use labeled routing. This result was improved to make routing table 
sizes smaller in [CGMZ05]. The problem of named routing over graphs with small doubling dimension has been 
studied in [KRX06] and [AGGM06], and references therein. To the best of our knowledge, there has been no prior 
work on dynamic graphs over doubling metric spaces and on control traffic overhead. It is worth pointing out that 
there is no direct correspondence between control traffic and memory. Bounds on memory do not take into account 
the amount of information which needs to be sent around in the network in order to build routing tables. A good 
illustration is the computation of the shortest path between two nodes u and v in a graph. While it is sufficient for 
every node on the path between these two nodes to have one entry for v (of roughly logn bits i.e., the name of 
the next hop), computing that shortest path requires a breadth first search of the communication graph and leads 
to a control traffic overhead of 0(n logn) bits. 

II. Models and Definitions 

A wireless network consists of a set of n nodes spread across a geographic area in the two-dimensional plane. 
We model the network region as the square area [0, \/n) x [0, y/n) . The n nodes move randomly in this area and 
we denote by (u) the position of node u at time t. The connectivity between two nodes is represented by an 
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edge on the connectivity graph Qn if they can communicate directly over the wireless channel. The connectivity 
between two nodes depends on the distance between the two nodes (and could also depend on the presence of other 
nodes, see Section III-Bb . We consider that when a node u transmits on the wireless channel, it broadcasts to all 
its neighbors in the connectivity graph Q$ . Consequently, one transmission of a packet is sufficient for all direct 
neighbors to receive that packet. To make the notation lighter, we will only add the dependence on time if it is 
necessary to avoid confusion. The distance d^'(u,v) between nodes u and v is the shortest path distance between 
these nodes in QrJ ■ Note that .) is a metric on Qn\ i.e., the distance between a node and itself is zero, the 
distance function is symmetric and the triangle inequality applies. We will now define a ball of radius R around a 
node u. It is simply the set of nodes within distance R of u. More formally, we can define it more generally for 
any metric space as follows: 

Definition 1: A Ball B^(u) around node u at time t in a metric space X is the set \v E X\d^ (u, v) < /?}. 
In order to bound the control traffic overhead, we will recursively subdivide the connectivity graph into balls. It 
will be crucial for us to bound the number of balls of radius R necessary to cover a ball of radius 2R around some 
node u. In other words, we want to find the smallest number of nodes V{ such that all nodes within 2R of u are 
also within R of some node vi. The notion of doubling dimension of a metric space captures this idea. 

Definition 2: The doubling dimension of a metric space X is the smallest a such that any ball of radius 2R can 
be covered by at most a balls of radius R, for all R > mm^ u ^ d(u, v) i.e., Vu G X 3 S u C X, \S U \ < a and 

£ U 4%) 

jes u 

Moreover, if a is a constant, we have the following definition: 

Definition 3 (Doubling metric space): We a metric space X is doubling if its doubling dimension is a constant. 
A good way to illustrate and understand the concept of doubling dimension and doubling metric space is to look 
at the metric space defined by a set of points X in R 2 with the Euclidean distance. A ball of radius 2R around 
a point x will simply be a disc of radius 2R around this point. To cover this disc, we will select a set of points 
such that all the surface is covered by the corresponding set of discs of radius R. Note that the number of discs 
required will not depend on R, and consequently this metric space would be doubling (see Figure \Q. Further, a 
metric space is said to be doubling if its doubling dimension is a constant, independent of the number of nodes n. 




Fig. 1. The metric space defined by a set of points in R 2 and the Euclidean distance is doubling. Indeed, we can cover a disc of radius 
2R by a constant (8 in this case) number of discs of radius R, whatever the value of R. 

In Section III- A I we describe the geometric random graph model, which will be the canonical model we will use 
to illustrate the ideas of the paper. We also give an example of a non-homogeneous network to which our results 
can be applied. In Section III-BI we will develop the model where connectivity is determined by the SINR, and we 
have uniform transmit power and full connectivity. We give the requirements for the mobility model to result in a 
smooth sequence of wireless network graphs in Section III-CI We state the underlying assumptions and give a table 
of notations in Section III-DI 
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A. Geometric random graph 

We denote the geometric random graph by Q(n,r n ) and define its connectivity as follows. 

Definition 4: A random geometric graph Q(n,r n ) has an unweighted edge between nodes u and v if and only 
if | \x(u) — x(v)\\ < r n , where {x(u)} are chosen independently and uniformly in [0, y/ri) x [0, y/n) . 
In this paper we will be interested in fully connected geometric random graphs, and therefore focus on the case 
r n > y/Togn [GK98]. As a natural extension, we can also define a sequence of random graphs ^^(n, r n ) with 
an unweighted edges between u and v at time t if ^^(ti) — a;(*)(i>)|| < r n . Whether each graph in the sequence 
Q^(n,r n ) corresponds to a random geometric graph as in Definition |U depends on the mobility model for the 
nodes. We discuss this in more detail in Section Hl-CI 

In Figure |2j we illustrate a non-homogeneous random network where connectivity is not completely geometric 
as in Definition [4] An obstacle prevents communication between neighboring nodes, and therefore illustrates the 
complexities of wireless network connectivity. This example is revisited in Section JIIJ where we show that though 
this connectivity graph is more complicated than Q(n,r n ), it is still doubling, and therefore the algorithms developed 
in this paper are applicable. This illustrates the advantage of our approach to network modeling. 




Fig. 2. n nodes are distributed uniformly at random on a square area of side y/n. A wall of width r„/c is added, which only has a small 
hole in the middle. Again, we assume r n > y/\og n. Nodes cannot communicate through the wall. Finally, we remove the nodes below the 
wall, which leads to an inhomogeneous node distribution. 



B. SINR full connectivity 

Since the wireless channel is a shared medium, the transmissions between nodes interfere with each other. 
However, the signal strength decays as a function of the distance traveled, and therefore we can define the SINR 
for transmission from node u to v as, 



SINR 



P n \\x(u) - x(v)\\-P 



x\v) 



-3- 



(1) 



E w ^u,v P \\ X ( W ) 

where j3 is a distance loss (decay) parameter depending on the propagation environment, P n is the common transmit 
power of the nodes and Nq is the noise power. We can of course easily adapt this to have power control for the 
nodes. A transmission is successful if the SINR is above some constant threshold value q. For static nodes, just as 
in the case of geometric random graph, we assume that the node locations {x(u)} are chosen independently and 
uniformly in [0, y/n) x [0, y/n) . This model for wireless networks has been extensively studied in the literature (see 
[GK00], [KV02]). The authors base their analysis of the capacity of wireless networks on a TDMA scheme for the 
SINR connectivity model of (Q~|). We argue here that the structure of the resulting connectivity graph is identical to 
that of Q(n,r n ), for r n > ylog n. Therefore, the results we prove for Q(n,r n ), would also be applicable to such 
graphs. In practice, it is a non-trivial task to design a distributed scheduling protocol (MAC layer protocol) that 
mimics the behavior of this TDMA scheduler. However, these MAC layer implementation issues are far beyond the 
scope of this document (see for instance [MSZ06]). We only make the argument here that the connectivity graph 
resulting from such a TDMA scheme would yield the same behavior as a Q{n,r n ). 

We will subdivide the network into small squares of side s n = ' £ f. We need to show that if two nodes u and v 
are in neighboring small squares (and so have the guarantee that they can communicate under the Q(n,r n ) model 
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as we will see in the sequel), then there exists a TDMA scheme that allows them to communicate under the SINR 
connectivity model of £[]). If this is the case, then we can apply the same proof techniques for both models. We let 
the maximum transmission power grow in the same way as we did for the Q(n,r n ) modejl i.e. P n < (N Q <;r n )P. 
Additionally, we want to design a TDMA scheme such that the capacity of all links is at least O(j^p^) [bits /sec]. 
It can be shown (see [RS98]) that every small square contains at most O(logra) nodes. Hence, we ask that the 
traffic can flow at constant rate independent of n between neighboring small squares, and that each node is tteated 
equally. Note that this requirement is very similar to the scheme proposed in [GKOO] in which one node per small 
square can transmit at constant rate to any neighboring squared. 

Theorem 1: There exists a TDMA scheme such that all nodes can communicate with any node located in a 
neighboring small square at a rate of Q( log ( n ) ) [bits/sec]. Hence, the aggregate traffic can flow between neighboring 
small squares at a constant rate independent of n. 

Proof: We take a coordinate system, and label each square with two integer coordinates. Then we take an 
integer k, and consider the subset of squares whose two coordinates are a multiple of k (see Figure [3]). By translation, 
we can construct k 2 disjoint equivalent subsets. This allows us to build the following TDMA scheme: we define k 2 
time slots, during which all nodes from a particular subset are allowed to transmit for the same duration of O(j^) 
seconds. Each small square contains at least one and at most O(logn) nodes w.h.p. (see [RS98] and the proof of 
Theorem [3). We assume also that at most one node per square transmits at the same time, and that they all transmit 
with the same power P n . Let us consider one particular square. We suppose that the transmitter in this square 
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x: nodes transmitting in this timeslot 



Fig. 3. Illustration of the TDMA scheduling scheme. 



transmits towards a destination located in a square at distance at most 1. We compute the signal-to-interference 
ratio at the receiver. First, we choose the number of time slots k 2 as follows: k = 4. To find an upper bound 
to the interferences, we observe that with this choice, the transmitters in the 8 first closest squares are located 
at a distance at least 3 (in small squares) from the receiver (see left-hand side of Figure [3]). This means that the 
Euclidean distance between the receiver and the 8 closest interferers is at least 2s n . The 16 next closest squares are 
at distance at least 7 (in small squares), and the Euclidean distance between the receiver and the 16 next interferers 

2 Note that the Q{n, r n ) model corresponds to the SNIR model without interferences. Indeed, if we remove interferences, two nodes can 

P„\\x(u)-x(v)\\-l 3 



communicate whenever 



No 



> <; for some threshold value <;. Hence, two nodes can communicate whenever — < 

3 The throughput achieved by this scheme is Q( -^=j==) [bits/ second/ node] when n source destination pairs are chosen uniformly at 
random. 



( ■ In particular, we let P n 
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is therefore at least Qs n , and so on. The sum of the interferences I can be bounded as follows: 

I <ZZi8iPn[2s n (2i-l)]- f3 

= {N oq r n f [2^]- /3 E£i^[(2 i -l)]-' 3 

This term clearly converges if /3 > 2. Now we want to bound from below the strength of the signal received from 
the transmitter. We observe first that the distance between the transmitter and the receiver is at most y/2{s%) < 2s n . 
The strength 5 of the signal at the receiver can thus be bounded by 

5 > P n min{l,2 S - /3 } 
= 0(1) 

Finally, we obtain the following bound on the SINR: SINR > N S +I - As the above expression does not depend 
on n, the theorem is proven. ■ 

C. Uniform speed-limited (USL) mobility 

Nodes are mobile and move according to the uniform speed-limited (USL) model, a fairly general mobility 
model defined next. The USL model essentially embodies two conditions: (i) the node distribution at every time 
step is uniform over the network domain, and (ii) the distance a node can travel over a time step is bounded. We 
restrict ourselves to the case in which the maximum speed is not dependent on n. In practice, of course, such an 
assumption is realistic since the maximum speed of the nodes will not increase when new nodes join the network. 

Definition 5: A collection of n nodes satisfy the uniform speed-limited (USL) mobility model if the following 
two conditions are satisfied: 

(i) At every time t, the distribution of nodes over the network domain is uniform; 

(ii) For every node u and time t, the distance traveled in the next time step is bounded, i.e., | |x( <+1 ) (m)-iW (u) \ \ < 
S. 

The USL mobility model is quite general. For example, it includes the following cases: (i) The nodes perform 
independent random walks with bounded one-step displacement. The random walks can be biased, and the dis- 
placement distribution does not need to be homogeneous over the node population. We have to assume that the 
nodes operate in the stationary regime, (ii) The nodes follow the random waypoint model (RWP). The system has 
to be in the stationary regime, (iii) The generalized random direction models from [SMS06], which interpolate 
between the random walk and the random waypoint cases, through a control parameter that can be viewed as the 
"locality" of the mobility process, (iv) We can also allow for models where nodes do not move independently. As 
an illustrative example, assume we uniformly place nodes on the square; the nodes then move in lockstep according 
to any speed-limited mobility process, maintaining their relative positions to each other. Observe that the uniform 
distribution is maintained for all time steps (note that we move on a torus), and that the speed-limited property is 
true by definition. 

We see that the USL class of mobility models is fairly general, and includes many of the models that have been 
proposed in the literature. For simplicity, we consider that time is discrete. In other words, we look at a snapshot 
of the network every AT seconds. At every time step, the connectivity between nodes will be modified. Hence, 
we will work with a sequence of connectivity graphs. In order to design a routing algorithm with a low control 
traffic overhead, we will need to understand how fast the graph distances between nodes can evolve over time. In 
particular, consider two nodes u and v at distance d = d l (u, v) at time t. We want to bound the multiplicative 
factor by which this distance can change in k time steps. Formally, we define k(t, d) as follows: 

Definition 6: We say that a communication network is k(t, d) -smooth if the shortest path distance between any 
two nodes u an v at shortest path distance d cannot change by more than a factor k(t, d) in r time steps i.e., Vit, v, 
we have: 

max \ d(^) ( U , v y 3®m ;- kM) 

Additionally, we simply say that the network is n-smooth if there exists a constant v such that K(ud, d) < k(u) = k 
independently of d. In this case, the distances grow at the same speed at all scales. In the sequel, we will bound 
K and v for our model. This USL property holds for a general class of random trip mobility models studied in 
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[BV05], where it is shown that the stationary distribution of such mobility models is uniform and ergodic. We 
restate this theorem without proof. 

Theorem 2: ( [BV05]) The random-trip mobility model has uniform stationary distribution on [0, a) x [0,a). 

D. Assumptions 

We consider that a time step AT is much larger than the round trip time (RTT) through the network i.e., the time 
scale for mobility is much larger than the time scale for communications. For clarity and in order to simplify the 
analysis, we will make the assumption that nodes can communicate instantaneously through the network. We also 
make the assumption that there is a random permutation it on the nodes, and that all nodes in the network know 
their rank in the permutation. In Section [VT] we will then drop these assumptions and consider practical aspects of 
the implementation. Finally, we say that a result holds with high probability (w.h.p.) if it holds with probability at 
least (1 — O(^)), for some constant p > 0. In Table U we summarize the notations used in this paper. 



d w (u,v) 

r n 
G(n, r„) 

k(t, d) 



Position of node u at time t 
Shortest path distance from u to v at time t 
Wireless communication radius 
Random geometric graph 
Ball of radius R around u 



TABLE I 

Table of notations 



III. Network Properties 

In this section, we prove some properties of the network models presented in Section [TTJ which are necessary 
to analyze the performance of our algorithm. We focus our attention on the geometric random graph Q(n,r n ), 
but all the arguments can be extended to the SINR full connectivity model with TDMA scheduling, discussed in 
Section Hl-BI In particular, for Q(n,r n ), we now consider the case in which the communication radius r n is such 
that r n = y/ (1 + e) log n > log 1 / 2 n, where e > 0. 

For uniform speed-limited (USL) mobility models discussed in Section III-CI at each time, the node locations 
{jW(u)} have a distribution that is uniform over [0, y/n) x [0, \/n) . Therefore, we now discuss the property of a 
sequence of geometric random graphs, G^(n,r n ), under USL mobility model. We subdivide the network area on 
which the nodes live into smaller squares of side where c is a constant chosen such that nodes in neighboring 
squares are connected (see Fig. |4| and that an integer number of squares fit into the network area. We arbitrarily 



Fig. 4. Nodes in neighboring squares are connected 

set c = \/5. We number the small squares from 1 to m = ^ r " e ^ 2 = ( 1+ ") C i ogn an d denote by E{ the event that 
small square i does not contain any node, in a sequence of length n p time steps, for some constant p. In the next 
theorem, we show that when nodes move according to USL mobility model, all small squares will be populated 
w.h.p. 

Theorem 3: There exists a constant p such that if we divide the network into small square of side ^ (with 
r n > y/log n), at every time step in a sequence of length n p , every small square contains at least one node w.h.p. 
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Proof: Consider a sequence of length Z 
j. Let m = 4-. We can compute: 



n' 



. Denote by E- the event the small square i is empty at time 



V 



ULU" 



E 



U) 



<zyJIUv[e\ 

2 nc 2 ! 

= z , nc e 

(l+e) logn 

< 7 nc2 1 
— (l+£)lognn( 1 +«) 

= Z-, 



(l+e)n e log n 

We can now choose p such that e — p > and the result follows. ■ 
It is immediate that a in single instantiation of the connectivity graph (i.e., time step), every small square is populated 
w.h.p. 

Corollary 1: With probability at least (1 — O(^)), there is no empty small square in a sequence of length 1. 

We are now ready to show that at every time step in a sequence of n p connectivity graphs, the connectivity graph 
is doubling w.h.p. Since we have a USL mobility model, any graph Q^(n, r n ) is statistically identical to G(n,r n ). 

Theorem 4: Q(n, yj (1 + e) log n) are doubling w.h.p. 

Proof: By Lemma [T] all small squares contain at least one node w.h.p. Consequently, neighboring squares 
(vertically and horizontally) have at least one communication link. Denote by C{m,r) the grid having the small 
squares as vertices, and with edges between vertical and horizontal neighbors. Consider a ball B upper = J3^ R (u) 
centered around some node u. Clearly, all nodes in B upper must be contained in a square which is part of B% Rc (u) 
i.e., B upper C B^ Rc (u). This follows from the fact that no node in B upper can be further away from u than 2Rr in 
Euclidean distance, and that the grid is fully connected w.h.p. Similarly, one can see that B R (u) C Bi ower = B R (u). 
This is a consequence of the fact that £ is a subgraph of Q, i.e., two nodes in small squares R hops a part in C 
cannot be more than R hops apart in Q (see Fig. [5]). For an appropriately chosen constant a, we have: 



B upper C B% Rc (u) C (J B&vj) C(JB R 



R\ U JJ 



(2) 



3=1 



3=1 



and Q(n, y/ (1 + e) logn) is doubling. ■ 
Note that it is possible to build a deterministic geometric graph for which this property does not hold (see Appendix 
I.O.al) . Further, one can show that Q(n,r n ) are not doubling w.h.p when r n < \/log n. We prove this result in 
Appendix I.O.bl At this point, we would like to emphasize that even though we analyze networks in which the nodes 
are uniformly distributed on a square area, the doubling property is a much more powerful tool. Indeed, our results 
and algorithms depend only on the doubling constant. Consequently, the algorithms and the bounds can be applied 
to any other type of networks or node configuration which lead to a doubling connectivity graph. For instance, one 
can consider the network shown in Figure |2 described in Section ITl-AI It can easily be shown by using a technique 
similar to the one used in Theorem [4] that this network is doubling. While we can seamlessly apply our routing 
algorithm to such a network, any classical geographic routing algorithm would fail or require a high control traffic 
overhead to get out of dead-ends. This is because nodes would get stuck against the wall when routing packets 
from the lower to the upper part of the network. In turn, this would considerably degrade the performance in terms 
of stretch and control traffic overhead with respect to the same network without a wall. In the next subsection we 
prove a set of sufficient conditions for a wireless networks to have a constant doubling dimension. 



A. Inhomogeneous Topologies 

In the first part of this subsection, we show that under certain conditions, the presence of topological holes 
(obstacles) in the network does not increase the doubling property, or only by a constant factor. In particular, we 
are interested in how we can alter the topology of a fully connected and dense network by removing nodes while 
still preserving the doubling property. In the second part, we will generalize this idea to arbitrary metric spaces. 
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Fig. 5. Proof of theorem [4] 



Consider a Q(n,r n ) with r n > y/log n, such that full connectivity is guaranteed. The network area is divided 
into squarelets of side r -^, where c is chosen such that nodes in horizontally and vertically adjacent squarelets 
are guaranteed to be within communication range. We now arbitrarily select squarelets and remove all nodes they 
contain. We denote the new graph we obtain by Q n . We denote by C n the full grid where the squareletes are 
vertices and by TL n the corresponding grid in Q n i.e., the thinned out grid obtained by selecting only non-empty 
squarelets in Q n . In C n , we add an edge between horizontally and vertically adjacent squarelets (see Fig. [6]). In 
TL n , we first add a an edge between horizontally and vertically squarelets containing at least one node. Then, for 
every pair of squarelets containing nodes that can communicate directly, we add an edge of weight corresponding 
to the distance between those two squarelets in C n . We add the edges from the shortest to the longest one, and 
only if no path of the same length already exists in TL n . We can now define a topological hole as follows: 

Definition 7 (Topological Hole): A set of horizontally, vertically and horizontally adjacent empty squarelets in 
the graph TL n is called a hole if adding a (virtual) vertex in all of the squarlets in that set modifies the distance 
between at least two vertices in TL n . 

Let us denote by the k th hole (k=l,2,3,...). We define the perimeter p(V&) of Vk as 2 times the maximum 
distance between any two vertices on the border of the hole i.e, in squarelets adjacent to the empty squarelets 
defining the hole. Note that for all u, v, we have cP~^ n (s u , s v ) > dP n (u, v). 

Theorem 5: Let P = maxfep(Vjt). Then, the doubling dimension is upper bounded by 0(P 2 ). 

Proof: Consider a ball (u) centered at u in Q n . First, observe that B^r (u) C Box^^ c (u), where 

B° x 2Rc ( u ) l& tne b° x cen t ere d at the squarelet containing u in C n which contains all nodes at "maximum norm" 

2Rc (i.e., /oo-norm) from this squarelet. In other words, all nodes within 2R hops from u in Q n must be in a 

r 

squarelet contained in this box. We will now cover this box with smaller boxes Box n , r (s v ). We need 
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Fig. 6. Graphs Q n , C n and 7i. n . The network area is divided into squarelets of side ^f- such that nodes in horizontally and vertically 
adjacent squarelets are guaranteed to be within communication range. 
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Perimeter of 64 



Fig. 7. A network with holes 



16R jy 47 = 64c 2 7 2 such boxes at most. Consider the same small boxes in 7i n . Pick one non-empty squarelet 
s Vi in each such small boxes. Note that the maximum hop distances between two squarelets in such a small box 
in C n is at most ^. For each of these hops, we might have to make a detour of at most P steps. Consequently, 
the same two squarelets could be at distance at most 4p i n T~L n . Observe now that for any two nodes u and 

w contained in squarelets s u and s w respectively, we have cP~^ n (s u ,s w ) > dP n (u,w). For each squarelet s v ., 
we pick one node V{ contained in this squarelet. Hence, for all nodes w contained in this small box, we have 
dpn (vi,w) < d^ n (s Vi ,s w ) < —. By setting 7 = P, we obtain the claim. ■ 
We can extend this result to the case where the network can be divided into convex sets. We define a convex set 
in Q n with slack as follows: 

Definition 8: Let ^ be a set of nodes in Q n . Let 7i n ^ be the squarelets in 7i n containing at least one node 
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in ^. We say that the set \f is convex if Vit, t> G ^, <i nn (s^,^) = d n {s u ,s v ), where s u and s v are the 
squarelets containing u and t> i.e., there must be at least one shortest path inside the convex set. We say that the 
set ^ is convex with slack P if d n (s u , s v ) < Pd n (s u , s v ). 
We can now state the following theorem 

Theorem 6: Let (j, £2, Cq be a partition of the network into g convex sets with slack Pi, P2, P q respectively. 
Let peri,per2, ...,per q denote the perimeter of the convex sets. The doubling dimension is then upper bounded by 

Proof: In the proof of Theorem [5 
in a box Box^ n (s u ), where p = 2Rc. We can cover each convex set Q intersecting this box with at most 



we have shown that any ball of radius 2R around some node u is contained 



per, 



small boxes of radius R/APi, as shown in Theorem pit A slack of P implies that by selecting one 



node in each of the small boxes, all nodes in the convex set are within R hops of this node in Q n . If the box 
/» 

Boxp n (s u ) is partitioned into several convex sets, selecting 4( ^p- ) 2 nodes in each convex set Q intersecting 

p i 

this box will in turn guarantee that all nodes are covered. ■ 
In practice, this result implies that if we are given a decomposition of the network into convex sets, we can bound 
the overall doubling dimension given the doubling dimension of each set separately. Further, this result implies 
that networks that consist of a small number of convex areas, which can each contain arbitrarily many small holes, 
have a low complexity in terms of doubling dimension. We will now relate the "shapes" of a topological hole to 
the doubling dimension. In particular, we will show that one can relate the doubling dimension to the maximum 
number of connected components in any square subarea. 

Theorem 7: For any 7 > 2, the doubling dimension a is such that 

a < 4"/ 2 c 2 max | number of convex disconnected components with slack 7 in Box^J^ 
Box Rh ^ 

Proof: In the proof of Theorem [51 we have shown that any ball of radius 2R around some node u is contained 
in a box Box^ n (s u ), where p = 2Rc. In turn, we showed that by dividing this box into smaller boxes of side 
R/j, and by selecting one node in each box, we could cover the larger ball of radius 2R. Now, in each small box 
of side -R/7, the presence of holes might create several disconnected components. However, we know that inside 
each such component, we can cover any convex subset with slack 7 with one nodes. The result follows. ■ 
This last result gives us a characterization of the alterations we can make to a fully connected G(n,r n ) network, 
while only affecting the doubling dimension by a constant factor. In particular, we can remove nodes as long as we 
do not create too many convex and disconnected components in any square subarea. Note that we can still remove 
arbitrarily many nodes as long as we only create small holes. Theorems [5j [5J [5J imply that topologies such as the 
one shown in Fig. [8] have a constant doubling dimension. The results stated above are special cases of the more 
general result detailed in the sequel. Indeed, we can relate the doubling dimension in a metric space to the doubling 
dimension in another metric space if we know the distortion of the embedding that maps the points in one metric 
space to the points in the other metric space. The example above is a special case of that setup where we map the 
nodes of a graph to points in Euclidean space. Consider two metric spaces (X, d) and (X', d'), where d and d 1 are 
distance functions which define a metric on the sets of point X and X'. We could for instance consider the two 
metric spaces (X, ||.||) and (Tt, d(., .)) i.e., the points in the plane with the Euclidean distance and the nodes in the 
graph with the shortest path distance. A metric embedding is a bijective function <fi : X — > X' which associates to 
a point in one metric space a point in another metric space. 

Definition 9 (Distortion of an Embedding): A mapping <fi : X — > X' where (X, d) and (X, d') are metric spaces, 
is said to have distortion at most D, or to be a D-embedding, where D > 1, if there is a K G (0, 00) such that 
Vx,y G X, 

Kd(x,y) < <t(<Kx),<Ky)) < KDd(x,y) 

if X' is a normed space, we typically require K = 1 or K = -p. An embedding has distortion D with slack e if all 
but an e fraction of node pairs have distortion D under <fi. Additionally, one can loosen this definition by allowing 



4 A convex area of perimeter q can always be included in a square area of side q. 
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Fig. 8. A network with topological holes and a constant doubling dimension. The size of the large holes grow with n, but the network can 
be divided into a constant number of areas, each being convex with slack O(l) i.e., each of the convex areas contains only obstacles with a 
constant perimeter or that can only increase the distance between nodes by a constant factor. Note that even though the doubling dimension 
is low, greedy geographic forwarding of packets would fail as packets would get stuck in dead-ends against the holes. Squarelets containing 
no nodes are hatched. 



slack. The slack is said to be uniform if each node has distortion at most D to a 1 — e fraction of the other nodes. 
Finally, an embedding with distortion D and slack e is coarse if for every node u the distortion is bounded to a 
node a distance greater than r e = inf {rs.t \B*(u)\ > en}. 

The doubling dimension of a metric space embedded into another metric space can be bounded as follows: 

Theorem 8 (Bounding the Doubling Dimension): Consider a metric space (H,d) embedded in another metric 
space (£, d!) by a function (f>. Let the doubling dimension of £ be /?. Let the distortion of this embedding be D. 
Then, H has doubling dimension a with a < 0((2L>) log/3 ). 

Proof: Choose any node u € 7i. If the above condition is fulfilled, the images of all nodes in B~2 R (u) can 
be at distance d 1 at most 2KDR from u at <j)(u). Hence, ^(B^iu)) C B2 KDR {4>{u)). We will now try to cover 
(p(B^ R (u)) by as few balls B RK (4>(v)) as possible (see Fig. |£l which illustrates this setup in the case when T~i is a 
graph and £ the Euclidean space). To do so, let us cover B2 KDR ((j)(u)) by small balls of radius KR in £. Covering 
B2 KDR (4>(u)) will require at most (3 log2D balls of radius KR in £, given that £ has doubling dimension j3. We 
know that d(u, v) < d'(u,v)/K, by definition |9l Consequently, <fi~ 1 (B RK (4>(v))) C B R (v). We can conclude that 

^c«) cupr^). ' ■ 

The presence of large obstacles in the network does not necessarily imply that the network is not doubling. In 
particular, 

Theorem 9: Consider a metric space £ with doubling dimension j3. A metric space TL that can be divided in 
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Fig. 9. Proof of Theorem [8] 



k sets Si, S2, —, Si-, such that each set embeds individually with distortion Di into £ has doubling dimension at 
most Z k J= iP 2log2Dj - 

Proof: Consider any ball of radius 2R in TL, such that the nodes in the ball belong to at least two different 
sets (otherwise the theorem is clearly true). Note that the radius of each of these subsets can be at most 4R. 
Consequently, we now that the part of the ball that belongs to Si can be covered by at most f3 2lo & 2D i (by applying 
Theorem [8] to cover a ball of radius AR by balls of radius R). The theorem follows. ■ 
We can now broaden the class of communication networks that have low doubling dimension. In particular, if 
we can subdivide the communication graph into a constant number of subsets, such that each one embeds with 
constant distortion into the Euclidean plane, the whole network is doubling. Consequently, topologies such as 
the one shown in Figure [10] are doubling. In this example, we embed an unweighted graph into the Euclidean 
plane. Note that the minimal Euclidean distance between nodes should be pr n (for some constant p), such that 
pr n d(u,v) < \\x(u) — x(v)\\ < 0(l)pr n . If this equation is true for all pairs of nodes, then the distortion is 
0(1). There is an issue when the nodes are neighbors in the communication graph, as the above rule implies that 
the Euclidean distance between such pairs of nodes should then be at least 0(r n ). However, we can ignore the 
distances below 2 as we will not cover balls of radius 1 (since we have a broadcast medium, the degree of a node 
does not impact the communication overhead). In such cases, it is obvious that geographic routing would fail, even 
though the inherent complexity of the network is low. Indeed, packets would get stuck against walls. Remarkably, 
our routing algorithm is oblivious to the topology and only depends on the doubling dimension. Hence, there is 
absolutely no need to detect or identify obstacles. The communication overhead will simply depend on the doubling 
dimension. 

B. Sequences of Communication Graphs 

In this subsection we study the behavior of a sequence of communication graphs, without any obstacles. We 
show that a sequence of G {t) {n,r n ) of length n p , for some constant p, with the USL mobility model is K-smooth. 
As already seen in Theorem [4] such a sequence of graphs is doubling at every time instant. 

Theorem 10: A sequence of QW (n, r n ) of length < n p , where nodes move according to the USL mobility model 
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Fig. 10. A set of doubling network topologies. The network is dense, and made inhomogeneous by the walls, which do no allow transmissions 
to go through. Note that the walls stretch when n grows, such that the network wide distortion also grows with n. Dashed lines indicate the 
separation into sets. 



with maximum constant speed S is 



max < 



J® ~ 7Z k „ ,V5V2(1 + - jm ) 



, VEV2 



smooth w.h.p. 

Proof: Consider two nodes u and v at Euclidean distance q£ = \\x u — x v \\2 at time t. Let qf^ = \ \x u — x v \\i = 
y^=i \ x m {u) — x m {v)\. Further, denote by d$ = d^\u, v) their shortest path distance at time t. One can see that 

< $t> < r q2 . Indeed, the shortest possible path will follow a straight line between u and v. The length 

of this line is and one hop can be of length at most r n . In the worst case, the shortest path from u to v will 
follow the shortest path in the grid formed by the small squares of side ^- = which exists w.h.p. Recall that 
we can only guarantee horizontal and vertical connectivity between small squares. The number of small squares in 



this path will be at most v ^ gl . One can easily show that qf 1 < y/2q^\ Let x = ( Xl ) = ( Xl j = (x u — x v ). 

q® = y/xj + xl = xiVl + s 2 



We have 



(l + s) Xl ^±f = q\ 



1+s 11 1+s 



Since, we have %y = ^-rjj- , the term is maximized when s = 1. In Figure [TT] we illustrate this point. Similarly, 

at time t + r, the shortest path distance will be bounded by ^ < dS t+T > < 92 . However, we know that 

the Euclidean distance can change by at most 2tS in r time step^l Consequently, 



q { 2 ] -^S < d(t+T) < V5V2(q? + 2rS) 
We can now bound the multiplicative stretch as follows: Hence, 



max 



max 



max \ r - d(t) , ^5^2(1 + \ = «(t, d) 

I 1 5^- 2 v / 5v^rS ^ J 
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Fig. 11. Upper and lower bounds for the shortest path 



One can now observe that the time it takes to multiply the shortest path distance between two nodes at distance d 
is proportional to d. Note that the larger the communication radius r n , the smaller n. Hence, the distance grows at 
most linearly with time. In particular, we have: 

Corollary 2: There exist constants u and k defined in the proof such that a sequence of n p connectivity graphs, 
under the USL mobility model with maximum constant speed S, is K-smooth w.h.p. 
Proof: By theorem [lOl we know that the sequence is 



max < 



,V5V2(l + 2rf ^- 



'-^-2^3 r " d '" 



-smooth w.h.p. Note that both terms decrease as a function of the communication radius r n . Hence, we can set r n = 1 
without decreasing k(t, d). Similarly, both terms go down when the distance e£w goes up. We can therefore also 
set d^> = 1, which is the smallest possible distance in an unweighted graph. Consequently, if we set r = vd® = v, 
we can now write 

k(t, d) < max { —, -^-^ — , vW2(l + 2vSy/Ey/2) } 



VEV2 



which is constant for v constant. 



IV. Routing Algorithm 

We develop the routing algorithm and its performance analysis for a general class of dynamic networks which 
produce a sequence of doubling and smooth connectivity graphs. We have seen in Sections [TT] and [III] that this 
applies to a class of wireless connectivity models with USL mobility. For notational convenience we illustrate the 
ideas for a sequence Q^\n,r n ) geometric random graphs with USL mobility. 

We decompose a time step into two phases: a beaconing phase and a forwarding phase. In the former phase, a 
set of routes are established by letting all or a subset of nodes flood the network at geometrically decreasing radii 
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and nodes register with beacon nodes. In the latter phase, this subset of routes is then utilized by source nodes to 
efficiently search for the destination. Every node is equipped with a routing table as shown in Table ITT] We first 



Node identifier 


distance [hops] 


level 


next hop 











TABLE II 
Routing Table RT 



describe two procedures used in the beaconing and the routing phase. 

flood(R,level) procedure: When a node u initiates the flood(R, level) procedure, it broadcasts a flood packet 
as shown in Table ITTT1 to its direct neighbors in Q(n,r n ). The hop count field is initialized to and the content 
of the Level field is set to the level of the beacon. How the level of the beacon is determined will be specified in 
the sequel. All nodes can compute the maximum hop count given the value of the level field in the packet. The 



Pkt. Type 


Node Id. 


Hop Count 


Level 


Oil) bits 


O(logn) bits 


O(logn) bits 


0(log A) bits 



TABLE III 
Flood Packet 



neighbors which receive this packet, after increasing the hop count by 1, add an entry to their routing table for node 
u if no entry for the same node and level with lower or equal hop count is present in the RT. The next hop field is 
set to the identifier of the node from which the packet was received. The level field in the routing table is set to the 
level given in the packet. In turn, the nodes which got the packet from u broadcast this packet to their neighbors. 
The latter follow the same procedure and update the routing table if necessary. The packet is discarded when the 
hop count reaches the maximum hop count (which is a function of the level). Note that with this procedure, every 
node forwards the packet at most once and the distance added to the routing table is the shortest path distance 
in Q(n,r n ). This procedure also allows us to establish a reverse path from all nodes that get the packet back to 
u. Indeed, it suffices for all these nodes to store the identifier of the node from which they received the packed 
Further, for any node v, that reverse path a shortest path to u. 

probe(relay,destination) procedure: This procedure consists in sending a probe packet (see Table HvT ) to a 
"relay" node for which the source has an entry in its routing table. The relay node will set the success bit to 1 if it 
has an entry for the destination and otherwise. We will make sure that all nodes on the path between the source 
and the relay node have an entry for the relay node in their routing table. Additionally, nodes on the path add a 
temporary entry for the source in the routing table. They set the next hop field to the identifier of the node from 
which they received the packet and leave the level and distance field empty. Upon receiving the packet, the relay 



Pkt. Type 


Relay Id. 


Dest. Id. 


Success 


O(l) bits 


O(logji) bits 


O(logn) bits 


1 bit 



TABLE IV 
Probe Packet 



node can either answer to the source on the reverse path we just created if the answer is negative. Alternatively, it 
can take action as explained in the sequel if it has an entry for the destination. 

We now separately detail the beaconing and the routing algorithms underlying our routing protocol 

A. Beaconing Algorithm 

In this subsection, we start by describing the first time step, when nodes have not yet moved and no information 
has been exchanged. In a static network, the information exchanged in this first step would be sufficient to setup 

6 The packet for which they modified their routing table 
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a complete routing infrastructure. On the other hand, in a mobile environment, we would need to cope with the 
dynamic topology and constantly update the routing tables. How we deal with a dynamic environment is explained 
at the end of the subsection. Let the cover radius at level i, for i = 1, ...,log A (A being the diameter of the 
network), be defined as r,j = 2* and the flooding radius at level i be defined as /j = n{ri + \ + r{), where k is a 
constant chosen such that n{vd, d) < nd, V<± In order for the routing algorithm to work properly, it is crucial that 
beacons at level i are within the flooding radius of the beacons at level i + 1, if they have common nodes inside 
their cover radii (see Fig. [T2l ). This is why we define the flooding radius above. Note that if the network is static, 




Fig. 12. The flooding radius is chosen in such a way that beacons at level i hear the floods of beacons at level i + 1. In static networks, 
this is n + n+i. 

we can set k to 1. Else, the value of k depends on how mobile the nodes are (see Def. [6]). 

The idea of the algorithm is to build a hierarchical cover of the network i.e., we would like every node in the 
network to be within rj hops of a beacon node at every level i. We say that when a node is within r, of a beacon 
b at level i, it is a member of b's cluster at level i, but it can only be in one cluster at every level. To achieve 
this, we let the nodes flood in a random order which can change at every time step. Every node u is a beacon at 
a given level /3(u). The flooding radius, however, will depend on the highest level at which a node is not covered. 
Let us denote by h(u) the highest level at which node u is not covered. When node u's turn to flood comes, it will 
determine the value of h(u) set j3{u) = h{u) and call flood(f h r u \,h(u)). A node v which receives this flood will 
determine the lowest level at which it could be a member of it's cluster, say l(v). That is, it will determine the 
lowest value j for l{v) such that d(u,v) < 2 J . This distance d(u,v) is known since v just received a flood packet 
from u. It will then become a member of it's cluster for all levels above l(v) for which it has no membership yet 
and are below P(u). If a node becomes a member of it's cluster, it sends a membership packet (see Table [VT> back 
to it. In this way, u learns the identifier of all nodes in its cluster. Note that u also applies this procedure to itself, 
and consequently could be a beacon at level i but not at level j < i. 
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Pkt. Type 


Node Id. 


Beacon Id. 


Level 


O(l) bits 


O(logn) bits 


O(logn) bits 


C>(log A) bits 



TABLE V 

Membership Packet 



The control traffic will be dominated by the messages sent back by nodes to beacons when they become members 
of a cluster, so should be rare. Moreover, we do not want the distance between a node and its beacon to grow by 
more than a constant factor. Since we assume that the maximum speed of the node is constant, the higher the level 
of a beacon, the more time it will take for nodes to double their distance to this beacon. We want to elect new 
beacons and update memberships only for levels at which the distances could have been multiplied by a constant 
factor. Recall that the network is K-constrained. Consequently, the distance S^> (u, v) between two nodes u and v 
cannot change by a factor k in less than vd time steps (see Corollary [2]). In particular, if a node is at distance 2 l 
of a beacon at the time it becomes a member of its cluster, then we have d t+u2 ' < k2\ Hence, we update the 
memberships at level i only every v2 l time steps (see Figure [T3l) . This will lead to a routing scheme in which the 



Layer 




d ■ 2 4 6 8... 2 J A 

Fig. 13. The memberships up to level i are updated every v2 r time steps. At the levels above, beacons elected at earlier time steps simply 
flood again. 

distances can be distorted by at most a constant factor to be calculated in the sequel. Additionally, in a dynamic 
environment, routes can break. This is why we let the beacons at all levels flood at every time step. Levels at 
which no membership updates take place simply use the floods of the beacons to update their routes toward theses 
beacons. This will ensure that a route always exists for all pair of nodes. 

In Figure [TJ] we give a simple with three levels. The beaconing algorithm is presented in Algorithm [TJ It is 
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Fig. 14. The example start with empty routing tables. First, on the left, node u\ floods at level 3. We focus on nodes ti2 and U3. Node u-i 
is within 8 hops from u\ but further away than 4 hops. Consequently, it can only had an entry for node u\ at level 3. At the same time, 
node uz can add an entry for node u\ at the levels 2 and 3, since it is at distance 4 of u\. Next, on the right, M2's turn to flood comes 
(right after Mi's turn). This node is already covered at level 3. Consequently, it will flood at level 2. The node U3 could potentially add an 
entry for this node at levels 1 and 2. However, it is already covered at level 2 and so adds only an entry for level 1. We do not show the 
entries beacons add for themselves. 



important to note that the routes are updated at every time step and consequently routing toward a beacon will 
always be successful. Further, when the membership at a given level i is updated, all the memberships at the levels 
j < i will also be updated, and all memberships at these levels canceled. 

B. Forwarding Algorithm 

The forwarding algorithms works as follows: a source node u with a message for a target node v searches for 
v by first probing all the level 1 beacon it knows of. To do so, it looks at its routing table and selects all nodes 
it knows of at level 1. If all answers are negative, node u probes all level 2 beacons it knows of. The procedure 
is repeated as long as all beacons answer negatively. A beacon at level i with an entry for the destination in its 
routing table does not answer directly to the source. Rather, this node will search downwards in the hierarchy by 
probing all the level i — 1 beacons it knows of. We show in the next section that one of these beacons must have 
an entry for the destination. That beacon in turn probes all the beacons in knows of at level i — 2. Meanwhile, the 
other beacons at level i — 1 will answer negatively to the beacon at level i. The procedure is repeated recursively 
until the target itself is reached. The target will then answer to the source on the reverse path which will later be 
used for communication between the source and the destination. To summarize, the forwarding algorithm starts 
with an "upstream" phase during which the source node probes beacons level by level until a beacon is found 
which has the destination in its cluster. That beacon then starte a "downstream" phase, during which we go down 
in the hierarchy. We illustrate the forwarding procedure conceptually in Figure [15] 
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Algorithm 1: Beaconing Algorithm at node u 
Data: Routing Table, Time t 
begin 

Let T = max {0 < j < log A\t modv2^ = 0}; 
Clear routing table entries with level < T; 
Let f](u) be the level at which u is a beacon; 
if ir(£) = u then 
if P(u) < T then 

Let h(u) be the highest level at which u is not covered; 

P(u) = h{u); 

end 

flood(f h{u) ,h(u)); 

end 

end 



C. Load-balancing 

This approach above guarantees a low network wide control traffic overhead. Even though over a long period of 
time all nodes will get approximately the same average overhead, beacons at the highest levels might get overloaded 
by the membership packets of the nodes in their cluster when a membership update takes place. These nodes will 
be hot spots in the network for a short period of time. To work around this problem, memberships can be distributed 
in the cluster instead of stored at the beacon itself. First, we now set f[ = K(2rj + i + r^). Additionally, whenever 
a beacon floods at level i, it includes its membership at level i + 1 in the packet. This information is stored by 
all nodes that receive this flood packet. This will guarantee that all nodes that are members of a cluster at level i, 
know how to reach all beacons at level i — 1 inside that cluster. A node that becomes a member of the cluster of 
beacon bi{u) at level i will now send its membership packet directly toward the beacon ^j-i(^) at level i — 1 inside 
this cluster with the identifier closest to n's. In turn, as soon as the packet reaches a node which is a member of 
^i_i(n)'s cluster at level i — 1, the membership packet is redirected toward the beacon ?/^_2(u) which is a member 
of ^i-i(u)'s cluster at level i — 1 and has the identifier closest to node u's. The process is repeated until we reach 
a single node, which will store u's identifier on behalf of bi(u). Note that the membership can only be registered 
at a single location in the cluster reachable through a unique sequence of clusters. This remains true even when 
nodes move. Indeed, the nodes in the cluster of bi(u) will only forward the packet to beacons at level i — 1 which 
were in the same cluster at the time the membership for this level got updated. Of course, whenever level j < i 
is updated, we do now not only need to send n's identifier toward its new beacon at that level. Additionally, the 
node that holds it's identifier at level j might not be reachable anymore through a path of clusters with identifiers 
closest to n's. Consequently, this node will need to forward n's identifier toward the beacon at level j with the 
identifier closest to n's. Again the process will be repeated recursively until a single node is reached. As we will 
see, the cost of avoiding hot spots is a factor log n in the total control traffic. Finally and most importantly, with 
this procedure beacons no longer get overloaded. Rather, the traffic is be distributed in its cluster. 

The data forwarding process remains the same except that the source node will not probe the beacon itself, but 
rather search for the node in the beacon's cluster that should hold the destination's identifier. If this node holds the 
identifier, it will then probe the beacons one level below in the same way. Recall that the nodes which potentially 
hold n's membership can be reached at any given instant in time through a unique sequence of clusters. The 
procedure is repeated until the destination is reached. 

V. Performance Analysis 

In this section, we analyze the performance of our algorithm analytically both in terms of control traffic and of 
route stretch. As in Section |IVJ we do this for a sequence of doubling and smooth connectivity graphs, and will 
use Q^\n,r n ) with USL mobility for illustration. 

The bounds derived in this section hold w.h.p. when we are in a sequence of length n p of a-doubling connectivity 
graphs. In the sequel, a, k and v are the constants derived in Section [Till Let us denote by A = 0(vT i<Wn.) )) 
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Fig. 15. Node u has a packet for node v. It searches in its routing table for all beacons it knows of at level 1 and sends them a probe 
packet containing v's identifier. None of the beacons at level 1 has an entry for this node and consequently they all answer negatively to 
node u. Next, node u repeats the same procedure with all the beacons it knows of at level 2. Again, all beacons answer negatively. On the 
third level, now, a beacon has an entry for node v. This beacon will probe all the beacons it knows of at level 2, while the other beacons at 
level three will answer negatively to u. A beacon at level 2 must have an entry for v. This beacon again probes all the beacons it knows of 
at level 1 among which one must have an entry for v itself. Meanwhile, the other beacons reply negatively as they do not have any entry 
for v. 



the diameter of the network. To bound the control traffic necessary for beaconing, we will rely on the a-doubling 
property of the metric space to show that a node can only hear a constant number of beacons at every layer. We 
will first show that a ball of radius 2R around any node u can only contain at constant number of balls (clusters) 
of radius R, when we select the centers of the balls of radius R in an arbitrary order and ensure that two centers 
cannot be closer than R. We will later use this result to show that a node can hear at most a constant number of 
beacons at any given level. 

Theorem 11 (Random Cover): Let B^iu) be a ball of radius 2R centered at it in a graph metric (X,d) with 
doubling constant a. Then, there exist at most k < a 2 nodes Vi, (i = 1,2,.., k) such that B* R (u) C (J f B R {vi) 
and min(i t j)d(vi,Vk) > R- 

Proof: By definition of an a-doubling metric space, there must exist a cover of a ball of radius 2R consisting 
of at most a balls of radius R. Recursively, there must also exist an -cover consisting of a 2 points. One can 
select at most one center Vi in each ball of radius R/2, as any other point inside this ball is within R of Hence, 
one can select at most a 2 such centers. ■ 
Corollary 3: Let B be a ball of radius R > R' in an a-doubling metric space (X,d). Then, one can select 
at most k < (jj?) 2l °9(<x) nodes v u (i = 1,2,.., A;) such that Bf (u) C U* s f'( u i) and min^divi, Vj) > R' . In 
particular, if R = r/R' for some constant R, then k is at most a constant {j]) 2l ° 9 ^ independent of n. 

Proof: Let R = 2 l R'. Hence, R' is doubled log^r times to obtain R. By Theorem [TT] B can be covered by 



23 



a 2log£ = (^)Mo»(a) baUs of radius R > a 

Here, one can think of the radius R of the large balls as the flooding radius, and of the radius R' of the small balls 
as the cover radius. Indeed, we use this result to show that a node u can hear the floods of all beacons within a 
given radius R. Moreover, this ball of radius R can contain at most (^ T ) 2lo 9i a ) beacons, since beacons must be at 
least R' apart. 



A. Control Traffic 

Theorem 12: The average control traffic overhead per time step for beaconing is at most 0(ralog 2 n) bits. 
Proof: We will analyze the control traffic at level i. Recall that a beacon at level i floods a distance /, = 
fv(2 4+1 + 2 % ) at every time step. Further, at the time the memberships are updated at level i, a beacon node at this 
level cannot be within rj = 2 l of another beacon at that level. If it were the case, this node would not elect itself as 
a beacon at this level. Level i is updated every u2 l time steps. Consider a node u. By Corollary [lOl no nodes that 
are further away than nfi hops at the time the memberships are updated at level i could move within /j of u in 
less than v2 l time steps. However, that is before this level is updated again. Consequently, the number of beacons 
whose flood can reach u at any given time step is at most the number of level i beacons in a ball of radius nfi at 
the time the membership is updated. In turn, node u will broadcast the flood packets of at most that many beacons 
for this level i. By Corollary d this number is a constant given by (^) 2lo 9(») = (3 K 2 ) 21oga . Given that there are 
O(logn) levels, that there are n nodes and that a flood packet has size Oilogn) bits, the average control traffic 
overhead per time step for beaconing is at most 0(n log 2 n) bits. ■ 
We now compute the control traffic overhead necessary for nodes to update their memberships with beacons. Recall 
that level i and all levels below are updated every v2 % times steps and that a node can only be a member of one 
cluster at every level. Furthermore, a node only becomes a member of a cluster if it is within 2* of the corresponding 
beacon. 

Theorem 13 (Membership Update Overhead): The average control traffic overhead per time step to update mem- 
berships without load-balancing is at most 

nlogAlogn 2 

= (J(n log n) 

v 

bits. 

Proof: Consider a sequence of T time steps. The memberships will be updated up to level i every v2 l time 
steps, so times in a sequence of length T. At the time of the update, a node can be at distance at most 2 l from 
a beacon at level i. Consequently, the overhead in bits generated by a node in a sequence of T time steps is upper 
bounded by Ya=i A ^2* log n = l ^sA i og n . ■ 
Finally, we will show that the average control traffic overhead when load-balancing is used is increased by at most 
a factor log n. 

Theorem 14 (Membership Update Overhead): The average control traffic overhead per time step to update mem- 
berships with load-balancing is at most 

ralog 2 Alog?i , 3 . 

= (Jin log n) 

v 

bits. 

Proof: Consider a sequence of T time steps. The memberships will be updated up to level % every v2 % time 
steps, so times in a sequence of length T. At the time of the update, a node can be at distance at most 2 l+1 from 
a beacon at level i — 1 inside its cluster at level i. Similarly, a node can be at distance at most 2 % from a beacon at 
level i — 2 inside its cluster at level i — 1. In the load balanced scheme, we have to count the overhead to go down 
the hierarchy of beacons. For a beacon at level i, this is at most 2 l x 2. Consequently, the overhead in bits generated 
by a node in a sequence of T time steps is upper bounded by 4 Y^i=i" 1°6 n = log n. However, node u 

is a member of a cluster at all log A levels. Recall that the node that holds u's identifier must always be reachable 
through a path by choosing the beacon (cluster) with the identifier closest to u's. Hence, whenever level % gets 

7 recall that when a node broadcasts a packet it is received by all direct neighbors in the connectivity graph. Consequently, there is one 
packet transmission per beacon of which a flood packet is received. 



In the load-balanced scheme, this constant is (5 



2 \ 2 log a 
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updated, all log A nodes that hold it's identity must follow the same procedure as u itself. We conclude that the 
overhead is upper bounded by log A4^j^ log n bits. ■ 

B. Route Stretch 

In this section we will show that the route found with the forwarding algorithm is only a constant factor longer 
than the shortest path. Additionally we show that the destination location discovery takes a negligible fraction of 
a flow throughput. 

Theorem 15 (Routing Stretch): The worst case multiplicative routing stretch is 0(1). 

Proof: We first analyze the stretch without load balancing. Consider that we want to route from a node it to a 
node v, and that we had 2 k < d(u, v) < 2 k+1 , the last time level k was updated before the route search takes place. 
Let us denote by bi(v) the beacon to which node v had registered the last time level i < k was updated before the 
route search takes place. Clearly, we have d(u,b v (k)) < n(2 k+1 +2 k ), and d(bi(v), &j_i(i;)) < k(2 % + 2 l ~ l ). This is 
true since the membership of node v at level i must have been updated at most v2 % time steps before the routing takes 
place, and that at the time the time level i gets updated, we have d(bi(v), &j_i(V))) < d(bi(v), v) + d(v, by 
triangle inequality. Note that d(bi(v), bi-i(v)) < /j_x and that d(u, b v (k)) < fk- Hence, a route must exist between 
it and v and the length A t+T \u,v) of the route at time t is at most: 

r W («.«) <£tiA = «£ti(2 i+1 + 2') 

= 3k Eti 2i = < 6nd® (u, v) 

In the worst cast, nodes u and v have moved closer together (by a factor k) while the beacons have moved further 
apart. Indeed, we have d^ t+T > (it, v) < nd^> (u, v) for r < v2 k as our network is K-constrained. Note that if we 
waited longer that u2 k , memberships would be updated again at level k and we could find another beacon at 

distance 2 k at most from v at level k. Hence, the worst case stretch is r J { ; + X' v> < 6k 2 = 0(1). ■ 
Every node can only hear floods from a constant number (/x = (3K 2 ) 2loga , see Theorem O of beacons at every 
level. Recall that the source will first probe all beacons at level 1, then all beacons at level 2 and so on. The 
procedure is repeated up to level k at which the source it will send a packet to bk(v). Note that the distance from 
it to this beacon can be at most n2 k+1 + 2 k = fk and so it must hear its floods. In turn, when routing down the 
hierarchy, beacon bj(v) will probe at most a constant number ((3K 2 ) 21ogQ of beacons at level j — 1. Finally, the 
distance between a node it and a beacon at level i can be at most fi and a probe packet will traverse at most 2f 
packets when a beacon at level i is probed (back and forth). This means that for discovery of the location of the 
destination, we need a probe overhead of at most fiQKd(u, v) packet transmissions. Therefore, this is a negligible 
part of the throughput of a flow since it consumes roughly the equivalent of a few packet headers of a flow from 
source to destination. A similar statement can be made for the load-balanced case. 

VI. Implementation Issues 

In Section |Vj we have computed worst case bounds which may be conservative in terms of constants. In this 
section, we explore this aspect by looking at simulation results for the control traffic and for the stretch. Recall 
that we had computed that for each of the O(logn) levels, a node has to retransmit a packet of at most (3K 2 ) 2Ioga 
beacons. Even if we set the maximum speed as well as the parameter v to 1, this is still \/T0 + 20 and consequently 
the constant in the bound on the overhead at least as high as (3(vl0 + 20) 2 ) 2 « 2.5 • 10 6 ! In Figure [T6l however, 
we show that in practice this constant is approximately 30. This simulation was run with 50 up to 10000 nodes 
moving at a maximum speed of 1. One can observe that the experimental scaling behavior corresponds extremely 
well to the theoretical behavior. To stress this fact, we also plot 100 log n as a benchmark. Note that the overhead 
is expressed in number of packets rather than bits (a packet being of size O(logn)). 

Similarly, in Figure [T71 we show that for a network of 1000 nodes, the stretch is at most 1.5 for all node pairs. 
If we compute the maximum theoretical stretch, we can show that it is again considerably larger and hence a 
pessimistic bound. These small constants could make a practical implementation realistic. 

We have made a certain number of assumptions in our models, which we now clarify. In practice, the random 
permutations on the nodes, which determines the order in which the flooding occurs, could be implemented by 
using random timers; more precisely, by letting all nodes draw a random delay independently of each other every 
AT seconds. Obviously, the interval from which nodes draw this delay should be made sufficiently large so that we 
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Fig. 16. Average control traffic overhead per node in packets as a function of the network size. Nodes move at a speed of maximum speed 
of 1. The confidence interval is given by the 95% and 5% percentiles. The size of a packet is O(logn) bits. We also plot 100 log n to show 
that our analytical predictions match the simulation results. 



can avoid collisions. However, a level in the hierarchy will be rapidly covered, and in a practical implementation the 
covers at different levels could be built in parallel. Further, different parts of the network are independent except 
at the highest level, and we could exploit this spatial diversity to parallelize the beaconing process. Hence, we 
speculate that it is possible to reduce the length of the beaconing phase to a small constant times the maximum 
round-trip time. Note that one could apply the algorithms to underlying networks that are not doubling. In this 
case, we would not be able to give provable bounds on the control overhead and the stretch as we did for doubling 
networks. 

VII. Conclusions 

In this paper, we show that a large class of wireless network models belong to a larger class of networks, the 
doubling networks, in which efficient routing can be achieved. To design an efficient routing scheme, one can 
hierarchically decompose the network by relying on the doubling property to prove that the control traffic overhead 
and the stretch will remain low, even for dynamic doubling networks. This holds for a fairly broad class of uniform 
speed-limited (USL) mobility models. One advantage of the proposed routing algorithm is that it is robust, in that it 
works well in certain situations in which other existing algorithms cannot work well. This was illustrated in Section 
III- A I for an example network with obstacles. We believe that many more such examples can be created where the 
use of the doubling rather than geographic properties would be crucial. To the best of our knowledge, our results 
are the first provable bounds for routing quality and costs for dynamic wireless networks. These techniques might 
give us insight into algorithm design for more sophisticated wireless network models. 
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Appendix 

a) Unit Disc Graphs: Another common model used in studies on wireless networks are Unit Disk Graphs 
(UDG), which are the deterministic variants of the random geometric graphs. The randomness of the positions of 
the nodes is removed and they can be placed arbitrarily on a finite of infinite area. The channel model is completely 
deterministic as before and nodes are connected if their Euclidean distance is below a threshold distance r, called the 
communication radius. In mathematical terms, two nodes u and v with positions x(u),x(v) G [0, R] 2 are connected 
if and only if \\x(u) — x(v)\\ < r. We will now show that there exist UDG which are not a-doubling (see Section 
|II]for a definition of an a-doubling metric). 

Theorem 16: There exists an infinite UDG for which is no constant that upper bounds the doubling dimension 
i.e., UDG are not doubling. 

Proof: Consider the graph shown in Figure [T8J To show that this graph is not a-doubling, we must show that 




Fig. 18. An infinite UDG obtained by deleting all the nodes in every second column of a grid, except for the nodes on the the middle row. 
Consequently, "columns" are 2r apart. 

there exists no constant such that all balls of radius R can be covered a constant a number of balls of radius R/2, 
for all R. Consider the ball centered around u in the figure. One can see that there are R/4 + 1 "columns" which 
cross the middle row at a distance less than R/2 from u (that is, the intersection of the column and the row is 
less than R/2 hops away from u). The intersection of each of these columns with B U {R) is of length more than 
R (see hatched zones on Figure [T8T). Consequently, for each of these columns there is at least one node at distance 
more than R/2 from the middle row. To cover these nodes, we need to place at least one ball of radius R/2 on 
each of these columns. Hence, the doubling dimension is lower bounded by R/4 and tends to infinity as R goes 
to infinity. ■ 
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One can notice that in the non-doubling UDG in the proof of Theorem [16] results from a careful construction. 
In Appendix I.O.bl we show however that such a structure will occur with high probability when r n < ylogn in 
random geometric graphs. 

b) Random Geometric Graphs with r n < \/log n: We first consider the case in which the communication 
radius r is such that r n = r = (log n) "2 < log 1 / 2 n and 9 G ] £, 1 ]. ( is a constant such that < £ < 1. 

Lemma 1: For any constant j3, there exists constants 7 > and b > such that a small square area of side 'yr 
with b nodes contains a subgraph of doubling dimension + 1 with probability p > 0. 

Proof: Consider the small square shown in Fig. [19] of side 77% where 7 is a constant independent of n to be 



Y 




Fig. 19. Graph for the proof of lemma Q] 

specified later. Subdivide the small square further into mini-squares of side r/c. Choose the constant c such that 
there exists a constant k satisfying y/2(k — 2) > c > \Jk 2 + 1. Under these conditions, two nodes in mini-squares 
separated by (k — 2) other mini-square will be connected, but not mini-squares r(k — 2)\/2/c apart (see right 
hand side of Fig. [T9l . Consider now the graph on the left hand side of Fig. [19] Assume that each full (colored) 
mini-square contains exactly one node. We now focus on the ball B2 R (u) and will lower bound the number of 
balls of radius R necessary to cover it. On the \_R/2\ first vertical branches from the left, the last node of the 
branch inside that ball (circled) must be covered by a ball of radius R centered on the same branch. This is clear 
since the length of the branch is larger than R. Consequently, the doubling dimension of this graph is at least 
L-R/2J > ^y^- We want the doubling dimension to be larger than f3, which can be easily achieved by choosing R 
such that > (3. Let R = 2/3 + 2 > 2/3 + 1. Further, we can now set 7 = (2R+5)(k - l)/c = (4/3 + 9)(fc- l)/c 

and b = (2R + 1) f 2 -^- + l] = (4/3 + 5) + 1 . This ensures that the doubling dimension is strictly larger 

than j3. 
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It remains to be shown that when such a small square contains b nodes, the graph constructed above occurs with 



probability p > 0. The number m of mini-squares contained in a small square of side is m = \t-^ = 7 c 

h T 

which is constant. Each node can fall in any of the m squares with equal probability. Hence, all m configurations 
are equiprobable and p = ^ > 0. ■ 
We number the small squares from 1 to m = -^p^ = 7 2 log 1 - 8 n an( ^ denote by -^i tne indicator variable that takes 
value 1 when small square i contains exactly b nodes. 

Lemma 2: There are at least n 1 / 2 squares containing b nodes with probability at least (1 — 0( ».^ B )) for n 
sufficiently large 
Proof: 

s[x] =£[ET=i^] 

E m ( n \( 1 Wl _ 1 \n—b 
i=l \b)\m) \ L m) 
> m (")6(J_)6(l_ ±)n 

> ^( T 2 log 1 " 9 n) h - 1 (1 - i)"^ 2 ^ » 

> ^(T'log 1 -^) 6 - 1 



> 0(n 1_0( ^ ) ) 

> 0(n s ) 

where S > | for n sufficiently large, since 9 > £. 

Let Si be the random variable representing the small square into which the i th node falls. Let F be the number of 
small squares containing exactly b nodes after all nodes have been placed. Then the sequence Zi = £ [F\S\, Si] 
is a Doob Martingale. One can show that F = f(Si,S2,--,S n ) satisfies the Lipschitz condition with bound 1. 
Indeed, changing the placement of the i th ball can only modify the value of F by at most 1. We therefore obtain: 

V[\F-£ [F]\ > n 5 / 8 ] < 2e- 2nl0/S " - " 



g 2n 1 / 4 



by the Azuma-Hoeffding inequality. Consequently, 



V [F < n 1 / 2 ] < V 



F< £\F]-n 5/8 



=n 7 /s_ n 5/s >n i/2 

< 2-3ttt < 2-4 



T7T 



and 



V [F > nV2] > (i _ 2 * ; 



e" 



It now remains to show that in this regime, Q(n,r) axe not doubling with high probability. 

Theorem 17: Q(n, (logn)2 _ 2), where 6 G ] £, 1 [ and £ is a constant such that < ( < 1, are not doubling with 
high probability. 

Proof: By Lemma [TJ for any constant (3, a small square area of side with b nodes contains a graph of 
doubling dimension > f3 with probability p > 0. By Lemma 12 there are n 1 / 2 such small squares containing b 
nodes w.h.p. Let F denote the number of small squares containing exactly b nodes. Consequently, the probability 
that at least one of this squares contains a graph of doubling dimension > [3 is given by: 

V [not doubling] = Y^=\ V [ not doubling |F = j]V [F = j] 

>(i-o(^))E™^(i-(i-ri J ) 

>(l_(l_p)»^)(l_ (_i„)) 

>(i-(i-p)« 1/2 ) 2 

>(l-^y) 

where x = (jz^) > 1. Consequently, with probability at least (1 — x0 2 (tl) ), there exists no constant which bounds 
the doubling dimension of Q(n, (log n)3 - f ). ■ 



